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METHOD AND SYSTEM FOR PREDICTING 
PHARMACOKINETIC PROPERTIES 

CROSS REFERENCE TO RELATED APPLICATION 
This application claims the benefit of U.S. Provisional Application No. 
60/21 1,864 filed June 14, 2000. 

Technical Field 

This invention relates to a method and system to predict pharmacokinetic 
(ADME) properties such as drug absorption (permeability), distribution, metabolism, 
and excretion, which are crucial properties in drug discovery. 

Background Art 

Experimental measurements to obtain pharmacokinetic properties are time- 
consuming and labor-intensive. Moreover experiments require a significant amount of 
actual compounds. Thus, the computational methods to predict such properties of 
virtual compounds are highly desirable in prioritization of targets prior to synthesis. 

So far, similar descriptors as conventionally employed in the quantitative 
structure activity relationship (QSAR) analysis (steric bulk, lipophilicity, HOMO 
energy, etc.) have been adopted in quantitative structure property relationship (QSPR) 
analysis to correlate with PK-related parameters (tl/2, clearance, or oxidation rate, etc.) 
(Lien, E. J. et al. Acta Pharm. JugosL 1984, 34, 123 — 131; Baeaernhielm, C. et al 
Chem.-Biol. Interact. 1986, 58, 277 — 288). Graph theory derived parameters 
(molecular connectivity indexes, etc.) have been also used for this purpose (Markin, R. 
S. et al. Pharm. Res. 1988, 5, 201 — 208; Garcia-March, F. J. et al. J. Pharm. 
Pharmacol. 1995, 47, 232 — 236). Recently reported QSPR methods necessitate 
calculation on 3D-structures that is still computationally intensive (Lombardo, F. et al. 
J. Med. Chem. 1996, 39, 4750-4755; Palm, K. et al J. Med. Chem. 1998, 41, 5382— 
5392; Clark, D. E. J. Pharm. Sci. 1999, 88, 815—821). The QSPR methods also 
necessitate complete set of molecular parameters (van de Waterbeemd, H. et al. Quant. 
Struct.- Act. Relat. 1996, 15, 480 — 490) that require experimental measurements to be 
determined. 
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2D-fingerprints are frequently employed in molecular similarity/diversity 
analysis (e.g. ISIS™/Base similarity search or SYBYL™/Selector), high-volume 
QSAR analysis (e.g. HQSAR, vide infra), and other drug discovery scenes. To date 
there has been no report on development of 2D-fingerprints descriptors to analyze 
pharmacokinetic properties. 

WO 98/07107 discloses a MOLECULAR HOLOGRAM QSAR (HQSAR™) to 
develop high volume QSAR models. HQSAR™ uses molecular hologram based on 
fragments counts to deal with mostly potency/activity. A symposium proceeding 
(Niwa, T. "Prediction of Human Intestinal Absorption of Drug Based on Neural 
Network Modeling"; 27 th Symposium on Structure- Activity Relationships held in 
Japan, Nov. 10, 1999) describes a method to estimate human intestinal absorption 
(HIA) based on molecular topological indexes derived from 2D-structure. 

It could be highly desirable to provide a system and method to predict 
pharmacokinetic properties of actual and virtual molecules with high 
performance (predictivity and speed) and wide applicability to diverse 
molecules. 

Brief Disclosure of the Invention 

This invention provides a new method and system for QSPR analysis and 
prediction based on only 2D-structure that allows us to predict hundreds of compounds 
rapidly. The method and system of this invention employs 2D-fingerprints, an array of 
the counts of functional groups as descriptors for QSPR. 

This invention provides a method for predicting pharmacokinetic properties of 
molecules comprising the steps of: 

(a) preparing 2D-structures of molecules used as a training set; 

(b) constructing a 2D-fingerprint by counting the number of structural descriptors 
that potentially relate to a pharmacokinetic property, either manually or automatically 
using internally developed macro; wherein said structural descriptors consist of 
predefined 20 to 80 atoms/fragments or substructures; 

(c) analyzing the obtained 2D-fingerprint by a statistical analysis method to correlate 
with the pharmacokinetic property of the molecule to yield a quantitative structure- 
property relationship (QSPR) model; and 



PC9980ACJG 



3 



(d) calculating the pharmacokinetic property of a trial molecule using the above 
obtained QSPR model. 

This invention also provides a system for predicting pharmacokinetic properties 
of molecules comprising: 
5 (a) means for preparing 2D-structures of molecules used as a training set; 

(b) means for constructing a 2D-fmgerprint by counting the number of structural 
descriptors that potentially relate to a pharmacokinetic property, wherein said structural 
descriptors consist of predefined 20 to 80 atoms/fragments or substructures; 

(c) means for analyzing the obtained 2D-fmgerprint by a statistical analysis method 
10 to correlate with the pharmacokinetic property of the molecule to yield a quantitative 

structure-property relationship (QSPR) model; and 
f** (d) means for calculating the pharmacokinetic property of a trial molecule using the 

if* above obtained QSPR model. 

SI Another aspect of this invention provides a method wherein the pharmacokinetic 

S j 1 5 property is absorption. 

fii 

J Another aspect of this invention provides a method wherein the pharmacokinetic 

property is distribution. 

pi Another aspect of this invention provides a method wherein the pharmacokinetic 

Q 

C t j property is metabolism 

*~ 20 Another aspect of this invention provides a method wherein the pharmacokinetic 

property is excretion. 

Another aspect of this invention provides a method wherein the internally 
developed macro comprises the macro script 2dfp.spl or 2dfp_abs.spl, written in 
SYBYL™ Programming Language (SPL). 

25 Preferably, each of the steps of the methods of the invention is carried out using 

molecular modeling software, databases or drawing software. More preferably one is 
such as SYBYL™, version 6.5 (Tripos Inc., St. Louis, MO). The database includes, for 
example, ISIS Base™ version 2.2.1 (MDL information Systems, Inc. San Leandro, 
CA.). The drawing software includes such as SYBYL™/SKETCH option, ISIS™ 

30 Draw version 2.2.1, Chem Draw Pro™ version 5.0 (CambridgeSoft. Corp. Cambridge, 
MA) and SMILES™ (Daylight Chemical Information Systems, Inc., Mission Viejo, 
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CA). Other modeling software, databases, and drawing software known to those of 
skill in the art can also be used. 

This invention enables us to perform virtual screening for synthetic targets and 
data mining using databases as well as drug design to optimize the pharmacokinetic 
5 profiles. Based on the QSPR model in this invention, it is possible to predict 
pharamacokinetic properties of molecules prior to synthesis, without labor-intensive 
and time-consuming experiment. This invention relies on 2D-fingerprint modeling 
requiring only 2D-structure, which enables us to perform rapid calculation to predict 
hundreds of compounds without tedious calculation about 3D-structure. Moreover, 
1 0 2D-fingerprint used in this invention comprises only 20-80 bits. 

Description of Figures 
p Figure 1 is a flowchart showing the overall process of the invention. 

1? Figure 2 shows a plot of actual vs. calculated log tl/2 . 

2! Figure 3 shows a plot of actual vs. calculated log(P £pp * 10 6 ). 

^1 1 5 Figure 4 shows a plot of actual vs. calculated logBB. 

Detailed Disclosure of the Invention 

O The term "molecules used as a training set" as used herein, refers to the 

p=j molecules whose pharmacokinetic properties have been already determined 

2i experimentally and used to develop a predictive QSPR model. 

h* 20 The term "pharmacokinetic properties" as used herein, refers to the properties 

of molecules related to metabolism, absorption (permeability), distribution, and 
excretion (ADME). 

A number of experimental methods or models are known in ADME. 
Examples of absorption studies are 1) kinetic studies based on measuring 
25 plasma concentration, urinary fecal excretions and gastrointestinal disposition after oral 
administration in vivo, 2) single-pass perfusion method, recirculation method, loop 
method in situ, and 3) everted sacs method, methods of using brush border membrane 
vesicles, isolated cells, and cultured cells (Caco-2) in vitro and the like. 

Examples of distribution studies are 1) the method of measuring concentration 
30 of target organs after administration by various technique such as HPLC, LC-MS, 
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autoradiography and microdialysis in vivo, 2) brain perfusion methods such as vascular 
reference method (brain uptake index) in situ, and 3) methods of using isolated cells or 
cultured cells (such as endothelial cell) in vitro and the like. 

Examples of metabolism studies are 1) kinetic studies based on measuring 
concentrations of drugs and the metabolites after adequate administration routes such 
as intravenous administration, administration per portal vain in vivo, and in situ, 2) 
kinetic studies such as a half-life of drugs in mammalian organ (li ver, kidney, intestine, 
etc. with slices, homogenates and microsomes etc) and in isolated cells or cultured 
cells such as hepatocytes in vitro and the like. 

Examples of excretion studies are 1) kinetic studies based on measuring 
concentration of drugs in urine, bile, feces etc after administration in vivo, 2) enzymatic 
studies of excretion via pumps such as P-glycoprotein, in vitro and the like. 

The term "2D-fingeiprint" as used herein, refers to a 2D-molecular measure in 
which a bit in a data string is set corresponding to atoms/fragments or substructures. 

The term "predefined atoms/fragments or substructures" as used herein, refers 
to atoms or functional groups relating to a phrmacokinetic property, which are based 
on the literature source (Bonse, V.G., Metzler, M. "Biotransformationen Organischer 
Fremdsubstanzen" (Yakubutu-Taisha) in Japanese Asakura, Tokyo (1980); Kato, R., 
Kamatani, T. "Yakubutu-Taishagaku" in Japanese, Tokyo-Kagaku-Dojin, Tokyo, 
chapter 4, 93-123 ,(1995)), otherwise refers to functional groups such as saturated or 
unsaturated bonds, rings (aromatic or cycloalkyl), amines, anilines, nittrogen in 
aromatics, imines/nitriles/guanidine/amidine, oxyamine(N-0)/nitro/azo-/hydrazin, 
amide/thioamide/sulfonamide/, alcohol/ether/aldehyde/ketone/ester/carboxylic 
acid/carbothioic acid/sulfinic acid/sulfonic acid, halogen, oxygen or sulfur functional 
groups, and total number of carbon, hydrogen, nitrogen, oxygen, sulfur or phosphorus 
atom. 

The term "internally developed macro" as used herein, refers to an internally 
developed Sybyl Programming Language (SPL) code. Preferable internally developed 
macro is as described in Working Examples 4 and 5. 

The QSPR model based on 2D-fingerprints for metabolism predicts half-life of 
molecules in a human liver microsome mixture with good predictivity. The 2D- 
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fingerprints for absorption are successfully employed to develop a higly predictive 
QSPR model on drug permeability across monolayer Caco-2 cells. Similarly the 
present 2D-fmgerprints/PLS modeling can be applied to develop statistically 
significant QSPR models on blood-brain barrier partitioning of the structurally diverse 
set. Thus, the method of this invention requiring only 2D-structures of the pertinent 
molecules enables to perform virtual screening of synthetic targets and data mining 
using molecular database as well as drug design to optimize the pharmacokinetic 
profiles. 

Figure 1 illustrates the method of this invention. This invention will be 
described in more detail with reference to Figure 1. Computational modeling studies 
can be carried out using molecular modeling software, preferably SYBYL™ on a 
Silicon Graphics Octane™ workstation. The method of this invention comprises the 
following steps: 

(a) 2D-structure of a molecule can be prepared by retrieving from a database 
such as ISIS™/Base, or by constructing manually with drawing software. The drawing 
software includes, for example, SYBYL™/SKETCH option (on the workstation), or 
ISIS™ Draw, Chem Draw™ and SMILES™ on (PC such as Windows NT client PC). 
The 2D-structure thus prepared can be transferred to the workstation, and stored in the 
molecular database. 

(b) The prepared 2D-structure of a molecule can be imported into molecular 
modeling software such as SYBYL™ as a MOL2 format. 2D-fingerprints can be 
constructed by the use of internally developed macro script 2dfp.spl or 2dfp_abs.spl, 
written in SYBYL™ Programming Language (SPL) implemented in SYBYL™, or by 
manually counting the number of the atoms/fragments or substructures. The macro 
program converts 2D-structures stored in the molecular database as a MOL2 format 
into a SYBYL™ line notation (SLN) format. Subsequently, the macro searches each 
SLN for the substructures potentially related to a pharmacokinetic property that match 
the queries described in the macro (as shown in Working Example 4), wherein the 
queries are predefined as the substructures (20 to 80 atoms/fragments). Finally the 
macro enumerates the substructure counts, and records them as 2D-fingerprints. 

(c) Statistical analysis is performed to obtain a correlation between the obtained 
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2D-fingerprints and the pharmacokinetic property. Any analytical method such as 
partial least square (PLS) algorithm, sample-distance partial least squares (SAMPLS; 
Bush, B. L. et al J. Computer-Aided Mol Design, 1993, 7, 587-619), genetic 
algorithm or neural network can be employed to yield an optimal quantitative structure 
property relationship (QSPR) model. 

(d) The pharmacokinetic property for trial molecules can be calculated based on 
the above obtained QSPR model. 

The pharmacokinetic properties of the molecule such as absorption, 
distribution, metabolism and excretion, can be apparent permeability coefficients 

(P app ) [cm/sec] blood-brain barrier partitioning ratio {(C brain /C blood ) = BB}, half- 

life(T!/2) in mammalian liver microsome and the like. 

The system of this invention can be constructed using appropriate computer 
hardware such as a Silicon Graphics Octane™ workstation and software as described 
above. 

This invention will be further described below with reference to the following 
Working Examples. 

Examples 
Example 1 

Development and validation of QSPR for half life in human liver microsome. 

Computational modeling studies were carried out using a Silicon Graphics 
Octane™ workstation. A congeneric series of 54 compounds of Formula (I)(as shown 
in the following Table 1 .) with a variety of substituent groups were used as a training 
set for analysis. 
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la cycloheptyl piperidinyl Ph 

2a cycloheptyl H 2 N(CH 2 ) 2 0- Ph 

3a cycloheptyl 4-aminopiperidyl Ph 

4a cycloheptyl H 2 N(CH 2 ) 2 C(0)- Ph 

5a cycloheptyl H 2 N(CH 2 ) 2 CONH- Ph 

6a cyclohepten-l-yl 4-aminopiperidyl Ph 

7a cyclooctyl H 2 NCH 2 CONH- Ph 

8a cycloheptyl H 2 N(CH 2 ) 3 - Ph 

9a cycloheptyl 4-aminocyclohexylamino Ph 

10a cyclohepten-l-yl piperazinyl Ph 

11a cycloheptyl piperazinyl Ph 

12a cycloheptyl H 2 N(CH ,) ? NH- Ph 
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# 


A Rl 


R2 






13a 


cycloheptyl 


H 2 NC(CH 3 ) 2 CH 2 NH- Ph 


5 


14a 


cycloheptyl 


N-methylpiperazin 


yl Ph 




15a 


cycloheptyl 


piperidinylamino 


Ph 




16a 


cycloheptyl 


4-aminopiperidyl 


CH 3 




17a 


cycloheptyl 


piperidinyl 


CH 3 




18a 


cycloheptyl 


H 2 N(CH 2 ) 10 NH- 


Ph 


10 


19a 


cycloheptyl 


4-aminoazetidinyl 


Ph 




20a 


cycloheptyl 


H 2 N(CH 2 ) 8 NH- 


Ph 




21a 


cycloheptyl 


(CH 3 ) 2 N(CH 2 ) 2 NH- 


Ph 




22a 


cyclooctyl 


N-methylpiperazinyl Ph 




23 a 


cycloheptyl 


piperazinyl 


isopropyl 


15 


24a 


cycloheptyl 


piperidinecarboximidamide Ph 




25a 


cycloheptyl 


H 2 N(CH 2 ) 6 NH- 


Ph 




26a 


cycloheptyl 


H 2 N(CH 2 ) 4 NH- 


Ph 




27a 


cyclononyl 


amino 


Ph 




28a 


cycloheptyl 


CH 3 NH(CH 2 ) 2 NH- 


Ph 


20 


29a 


cyclooctyl 


piperazinyl 


CH 3 




30a 


cycloheptyl 


4-aminopiperidyl 


vinyl 




31a 


cycloheptyl 


isopropyl 


Ph 




32a 


cycloheptyl 


2-guanidinoethyl 


Ph 




33a 


cycloheptyl 


mathanesulfoonyl 


Ph 


25 


34a 


cycloheptyl 


piperidinyloxy 


Ph 




35a 


cycloheptyl 


dimethylamino 


Ph 




36a 


cycloheptyl 


amino 


Ph 




37a 


cycloheptyl 


CH3CONH- 


Ph 




38a 


cycloheptyl 


hydroxypiperidinyl 


Ph 


30 


39a 


cycloheptyl 


H 2 N(CH 2 ) 3 S0 2 - 


Ph 




40a 


cycloheptyl 


methylamino 


Ph 




41a 


cycloheptyl 


methyl 


Ph 




42a 


cyclooctyl 


piperazinyl 


CH 3 




43a 


cycloheptyl 


isobutyl(NH 2 )CHCONH- Ph 


35 


44a 


cycloheptyl 


methylamino 


CH 3 




45a 


cycloheptyl 


methoxy 


Ph 




46a 


cyclooctyl 


methylamino 


normalprop 
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Table 1. (continued) 

# A Rl R2 

47a cyclooctyl methylamino CH 3 

48a cyclooctyl methylpiperazinyl CH 3 

5 49a cycloheptyl H Ph 

50a cyclononyl methylamino CH 3 

51a cyclononyl methylpiperazinyl CH 3 

52a cycloheptyl isobutyl(NH 2 )CHCONH- CH 3 

53a cycloheptyl H 2 N(CH 2 ) 2 CONH- Ph 

10 54a cycloheptyl H , N(CH ,), CCONH- Ph 

Half-life (tl/2) in vitro for each compound was determined by HPLC analysis 
of the reaction mixture with human liver microsome. The employed 2D-structures 
were retrieved from ISIS™/Base (version 2.2.1, MDL Information Systems, Inc., San 
15 Leandro, CA) or constructed with ISIS™/Draw (version 2.2.1, MDL Information 
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Systems, Inc., San Leandro, CA) on a WinNT client PC, followed by being transferred 
ffi to the Octane workstation and stored in a molecular database. The 2D-fingerprints 

were constructed by use of a newly developed macro script 2dfp.spl, written in 
SYBYL™ Programming Language (SPL), which was implemented in SYBYL™ 
20 (version 6.5, Tripos Inc., St. Louis, MO). The macro program converted 2D-structures 
stored in the molecular database as MOL or MOL2 format into SYBYL™ line notation 
(SLN) format, and counted the number of the atoms or functional groups that matched 
queries defined in a table described in the macro program. The atoms or functional 
groups susceptible to be involved in metabolism were assigned on the basis of the 
25 literature source (Bonse, V. G., Metzler, M. "Biotransformationen Organischer 
Fremdsubstanzen" (Yakubutu-Taisya, in Japanese) Asakura, Tokyo (1980); Kato, R.; 
Kamataki, T. "Yakubutu-Taisyagaku" in Japanese, Tokyo-Kagaku-Dojin Tokyo 
(1995)). Partial least square (PLS) algorithm in QSPR module in SYBYL™ was 
employed to correlate the aforementioned 2D-fingerprints and tl/2 to produce QSPR 
30 model. Thirty-eight bits out of whole 2D-fmgerprints used since 25 bits with all the 
same value or 0 were dropped. SAMPLS run in crossvalidation step (leave- 1 -out) 
identified the optimum PLS component as 5 (N = 54, Std. Error_prediction = 0.414; q 2 
= 0.518). Non-crossvalidation PLS analysis resulted in a significant five-component 
model with the following statistics: Std. ErrorJEst. = 0.219, r 2 = 0.865, F(nl = 5, n2 = 
35 48) -61.3. 

Figure 2 shows the plot of actual vs. calculated log tl/2 (closed circles). For 
validation of the present QSPR model, the prediction of half-life for the test set (12 
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compounds) was performed. As indicated open squares in Figure 2, the model has a 
fairly good predictivity, which allows us to prioritize the targets for synthesis. 

Example 2 

Development of QSPR for Caco-2 permeability. 

Unless otherwise noted similar computational molecular modeling were 
performed as described in Example 1. Table 2 enlists 21 structurally diverse 
compounds as a training set, whose apparent permeability coefficients (P app ) [cm/sec] 
of a compound across Caco-2 cells was used as in literature source (Yee, S. Pharm. 
Res. 1997, 14, 763 — 766). The counts of substructures to match with the predefined 
queries were encoded as a array of integers by a similar SPL script (2dfp_abs.spl) to 
afford 2D-flngerprints as descriptors employed in the correlation analysis. SAMPLS 
run in cross validation step (leave- 1 -out) identified the optimum PLS component as 2 
(N = 21, Std. Error_prediction = 0.444; q 2 = 0.463), Non-crossvalidation PLS analysis 
resulted in a significant two-component model with the following statistics: Std. 
Error_Est. = 0.254, r 2 = 0.824, F(nl = 2, n2 = 18) = 42.1. Figure 3 shows the plot of 
actual vs. calculated log(P app * 10 6 ). 



Table 2. Training set compounds with apparent permeability. 



Compd. 


P * 10 6 

r app 1 v 

(cm/sec) 


Compd. 


P * 10 6 

r app 1 w 

(cm/sec) 


Compd. 


P * 10 6 
1 app 1 v 

(cm/sec) 


Azithromycin 


1.04 


Diazepam 


70.97 


Prazosin 


43.60 


Benzylpenicillins 


1.96 


Erythromycin 


1.80 


Propranolol 


27.50 


Caffeine 


50.50 


Fluconazole 


29.80 


Quinidine 


20.40 


Chloramphenicol 


20.60 


Ibuprofen 


52.50 


Tenidap 


51.20 


Clonidine 


30.10 


Imipramine 


14.10 


Testosterone 


72.27 


Desipramine 


21.60 


Methotrexate 


1.20 


Trovafloxacin 


30.23 


Dexamethasone 


23.40 


Naloxone 


28.20 


Ziprasidone 


12.30 



Example 3 

Development of QSPR for blood-brain barrier partition. 

Unless otherwise noted, similar molecular modeling was performed as described 
in Example 1. Blood-brain barrier partitioning ratio, {log(C brain /C blood ) - logBB} for 
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"drug-like" compounds (N = 35, Chart 1) as a training set were used as in literature 
source (Lombardo, F. et aL, J. Med. Chem. 1996, 39 y 4750-^755.). The 2D-fmgerprints 
were calculated as above example. PLS modeling to correlate 2D-fingerprints with 
BBB partitioning ratio showed the following statistics. Crossvalidation (SAMPLS, 
5 leave-l-out): the optimum PLS component = 3, N = 35, Std. Error_prediction = 0.69; 
q 2 - 0.29. Non-crossvalidation: Std. ErrorJEst. = 0.38, r 2 = 0.78, F (3> 31) = 37. 4. 

Chart 1. Compounds employed in the analysis, (compound 36 for validation) 
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Chart 1 (continued). Compounds employed in the analysis, (compound 36 for 
validation) 
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Chart 1 (continued). Compounds employed in the analysis, (compound 36 for 
validation) 




OH 



Example 4 

SPL macro (2dfp.spl) to prepare 2D-fingerprints for half life. 

uims define macro 2dfp sybylbasic yes 

## 

## Set the Source Database, and Column-Names File. 

## 

setvar source %promptif("$l" "STRING" "MYFILE.MDB" "Source Database.mdb" 
"Database with molecules to be calculated") 

setvar resultsFP %promptif("$l" "STRING" "Columns.txt" "Filename storing column 
names" "Text file to store column names") 

## if%not(%mols(*)) 

## %dialog_message(ERROR "There are no molecules." "No Molecules") 

>$NULLDEV 
## return 
## endif 



# 



PC9980ACJG 




# Set the molecule area to calculate "2D-FingerPrint". 

# Note that Scurrent molarea is defined by the "calling" 

# table when adding a column of data. 

# 

localvar mol_area 
if $1 

setvar mol_area $1 
else 

setvar mol_area $current_molarea 
endif 

database open $ source read 

## 

## Loop over all molecules in the source database 

## 

for j IN %database(*) 

database get "$j" $mol_area 

# 

# set the SLN expression for the molecular area 
# 

setvar sln exp %sln($mol_area) 
setvar ARRAY 

# 

# # items + 1 (compd_num) BIT's will be used 
# 

############ @compdnum) ################## BIT 1 
setvar ARRAY[01] %mol_info($mol_area name) 

###### Exp-Generator_read(file_ID) looping is another choice 

################ Unsaturated bonds ########## BIT 2-4 
## fpla) Unsaturated bonds (aromatic) 
setvar query Any: Any 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 
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setvar ARRAY[02] $BIT 
## fp lb) Unsaturated bonds (bouble) 

setvar query Any=Any 
setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

setvar ARRAY[03] $BIT 
## fplc) Unsaturated bonds (triple) 

setvar query Any#Any 

setvar BIT %count(%search2D($ski_exp $query NoDuplicate 0 yes)) 
setvar ARRAY[04] $BIT 

################ ring (topology) ############# BIT 5-15 
## @fp2a) 3-membered ring 

setvar query Hev[l]~Hev~Hev@l 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 
setvar ARRAY[05] $BIT 
##@fp2b) 4-membered ring 

setvar query Hev[l]~Hev~Hev~Hev@l 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 
setvar ARRAY[06] $BIT 
##@fp2c) 5-membered ring 

setvar query Hev[l]~Hev~Hev~Hev~Hev@l 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

setvar ARRAY[07] $BIT 
##@fp2d) 6-membered ring 

setvar query Hev[l]~Hev~Hev~Hev~Hev~Hev@l 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

setvar ARRAY[08] $BIT 
##@fp2e) phenyl ring 

setvar query C[1]:C:C:C:C:C:@1 

setvar BIT %count(%search2D($sln __exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[09] $BIT 
##@fp2f) 7-membered ring 

setvar query Hev[l]~Hev~Hev~Hev~Hev~Hev~Hev@l 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

setvar ARRAY[ 1 0] $BIT 
##@fp2g) 8-membered ring 

setvar query Hev[l]~Hev~Hev~Hev~Hev~Hev~Hev-Hev@l 
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setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

setvar ARRAY[1 1] $BIT 
##@fp2h) 9-membered ring 

setvar query Hev[l]-Hev~Hev-Hev-Hev-Hev~Hev~Hev-Hev@l 
5 setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[ 1 2] $BIT 
##@fp2i) 10-membered ring 

setvar query Hev[l]~Hev-Hev~Hey~Hev~Hev-Hev~Hev~Hev-Hev@l 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
10 setvar ARRAY[13]$BIT 

##@fp2j) 1 1-membered ring 

setvar query Hev[l]-Hev~Hev-Hev~Hev~Hev~Hev-Hev-Hev~Hev~Hev@l 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[ 14] $BIT 



^ 15 ##@fp2k) 12-memberedring 

@g setvar query 

^ Hev[ 1 ]~Hev~Hev~Hev~Hev~Hev~Hev~Hev~Hev~Hev~Hev~Hev@ 1 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
01 setvar ARRA Y[ 1 5 ] SBIT 

y 20 



Q ################ Elements _Overall ############# BIT 1 6-22 

■JJ ##@fp3a) total Hetro atoms 

SJ setvar query Het 

p setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

^ 25 setvar ARRAY[16] SBIT ! 

##@fp3b) total Halogen 

setvar query Any[is=F,Br,Cl,I] 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[17] SBIT 
30 ##@fp3c) total N 

setvar query N 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[18] SBIT 
##@fp3d) total NH 
35 setvar query NH 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[ 19] SBIT 
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##@fp3e) total O 

setvar query O 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[20] $BIT 
##@fp3f) total OH 

setvar query OH 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[21] $BIT 
##@fp3g) total S 

setvar query S 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 
setvar ARRAY[22] $BIT 

################ Methyl, terminal ############# BIT 23-26 
##@fp4a) C-Methyl (omega-Oxidation) 
setvar query C-CH3 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[23] SBIT 
##@fp4b) N-Methyl (N-demethylation) 
setvar query N-CH3 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY [24] SBIT 
##@fp4c) O-Methyl (O-demethylation) 
setvar query 0-CH3 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[25] SBIT 
##@fp4d) S-Methyl (S-demethylation) 

setvar query CH3-S[F]-Any[NOT=H*] 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[26] SBIT 

################ Methylene -CH2- ############# BIT 27-30 

##@fp5a) Methylene group 

setvar query Any[NOT=H*,N,0]-CH2-Any[NOT=H*,N,0] 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY [2 7] SBIT 

##@fp5b) N-Methylene 
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setvar query N-CH2-Any[NOT=H*] 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[28] $BIT 
##@fp5c) O-Methylene 
5 setvar query 0-CH2-Any[NOT=H*] 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[29] $BIT 
##@fp5d) S-Methylene 

setvar query S[F]-CH2-Any[NOT=H*] 
10 setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

setvar ARRAY[30] $BIT 

################ Methine >CH-, Allylic/Benzylic H (to be absorbed) 
################ BIT 3 1-36 
15 ##@fp6a) Methine group 

setvar query Any[NOT=H* ) N,0,S]-CH(-Any[NOT-H*,N,0,S])- 

Any[NOT=H*,N,0,S] 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[31] SBIT 
20 ##@fp6b) Benzylic H (Ar-CH) (if Ph-CH2, then the count =2) 
setvar query CHC(:Any):Any 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[32] SBIT 
##@fp6c) Allyl H (if CR=CR-CH2, then the count =2) 
25 setvar query CHC(=C) 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[33] SBIT 
##@fp6d) N-Methine 

setvar query N-CH(-Any[NOT=H*])-Any[NOT=H*] 
30 setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[34] $BIT 
##@fp6e) O-Methine 

setvar query 0-CH(-Any[NOT=H*])-Any[NOT=H*] 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
35 setvar ARRAY[35] SBIT 

##@fp6f) S-Methine 

setvar query S-CH(-Any[NOT=H*])-Any[NOT=H*] 
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D 



Li 



setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[36] $BIT 

################ Nitrogen containing Compounds ############# BIT 37-49 
5 ################ Amines / Imines / Nitrile ############# BIT 37-46 
##@fp7a) Primary Amines, unbranched 
setvar query NH2CH2 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 
setvar ARRAY[37] $BIT 
1 0 ##@fp7b) Primary Amines, branched 

setvar query NH2CH(Any[NOT=H*])(Any[NOT==H*]) 
setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY [3 8] $BIT 
##@fp7c) Primary Amines, branched 
1 5 setvar query NH2C(Any|>rOT=H*])(Any[NOT=H*])(Any[NOT=H*]) 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY [3 9] SBIT 
##@fp7d) Primary Anilines (Ar-NH2) 
01 setvar query NH2C:Any(:Any[NOT=H*]) 

y 20 setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

□ setvar ARRAY [40] SBIT 

##@fp7e) Secondary Amines, 

setvar query NH(C[NOT=C=0])C[NOT=C=0] 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
25 setvar ARRAY[41] SBIT 

##@fp7f) Tertiary Amines 

setvar query N(C[NOT=C=0])(C[NOT=C=0])(C[NOT=C=0]) 
setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[42] SBIT 
30 ##@fp7g) Imines 

setvar query N=C 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[43] SBIT 
##@fp7h) Nitrile 
35 * setvar query C#N 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[44] SBIT 
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##@fp7i) N in aromatics 

setvar query Any[is=N,C]:N:Any[is=N,C] 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 
setvar ARRAY[45] $BIT 
5 ##@fp7j) Guanidine 

setvar query NC(=N)N 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 
setvar ARRAY[46] $BIT 

1 o ################ N-O, Nitro, N-N ####################### BIT 47-49 
##@fp7k) NO {Hydroxyamine, Oxime, Hydroxamic acid, ....) 
setvar query N-O 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

setvar ARRAY[47] $BIT 
~ 15 ##@fp71) Nitro (count =2), Nitroso (count =1) 

Qj setvar query N(=0) 

^ setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

Sj setvar ARRAY[48] $BIT 

f\ ##@fp7m) N-N 

20 setvar query N~N 

Q setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

2 setvar ARRAY [49] $BIT 

%** 

O ################ Amide, Ester, Sulfonamide ################ BIT 50-52 

fr ~' 25 ##@fp8a) Ester 

setvar query C(=0)OC 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[50] $BIT 
##@fp8b) Amide 
30 setvar query NC(=0) 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[51] $BIT 
##@fp8c) Sulfonamide 

setvar query NS(=0)(=0) 
35 setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 

setvar ARRAY[52] SBIT 
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################ Ketone, Aldehyde, Alcohol, Thiol, Sulfide ### BIT 53-59 
##@fp9a) Primary Alcohol 

setvar query CH2(OH)(~Hev) 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 
setvar ARRAY[53] $BIT 
##@fp9b) Secondary Alcohol 

setvar query CH(OH)(~Hev)(~Hev) 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

setvar ARRAY[54] $BIT 
##@fp9c) Ketone, Aldehyde 

setvar query Any[is=H,C]CC(=0)(Any[is=H,C]) 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 

setvar ARRAY[55] $BIT 
##@fp9d) COOH 

setvar query Any[is=H,C]CC(-0)(OH) 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 
setvar ARRAY [5 6] $BIT 
##@fp9e) Sulfide 

setvar query CS[F]C 

setvar BIT %count(%search2D($sln_exp $query NoDuplicate 0 yes)) 
setvar ARRAY[57] $BIT 
##@fp9f) Thiol 

setvar query S[F]H(C) 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[58] $BIT 
##@fp9g) Thiocarbonyl 
setvar query C=S 

setvar BIT %count(%search2D($sln_exp Squery NoDuplicate 0 yes)) 
setvar ARRAY[59] $BIT 

echo $ ARRAY 
echo $ ARRAY » SresultsFP 
zap Smol area 



endfor 

database close 
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## 

## Announce of the completion & the location of the results 
echo 

echo "job completed on %system(date)" 
5 echo "The results is stored SresultsFP as a text file," 
echo "please import if from MSS (table)." 

echo " => custom format, space separated, column label used" 



10 



15 



yy 

0 s 
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Example 5 

By use of a similar method described in example 4, 2D-Fingerprints for Caco-2 
permeability and blood-brain barrier partition were prepared based on the following 
description. 



Fp-ID 


Name 


Query 


Alkyl Amines 


fpla 


Primary 


NH2C[NOT=C=0 3 C=S ) C:Any,Any[IS=C,N]C(=N)N J C[l](Any[IS=0,S,N] 
Any=:AnyAny=@l),C[l](=AnyAny[IS=0,S J N]Any=:Any@l)] 


fplb 


Secondary 


NH(C[NOT=C=0,C=S,C:Any,Any[IS=C,N]C(=N)N,C[l](Any[IS=0 J S,N]A 
ny=:AnyAny=@l),C[l](=AnyAny[IS=0,S,N]Any=:Any@l)])(C[NOT=C= 
0,C=S,C:Any,Any[IS=C,N]C(=N)N J C[l](Any[IS=0,S,N]Any=:AnyAny= 
@l),Cril(=AnyAnyriS=0,S,NlAny=:Any@l)l) 


fplc 


Tertiary 


N(C[NOT=C=0,C=S,C:Any 5 Any[IS=C,N]C(=N)N ) C[l](Any[IS=0,S,N]An 
y=: AnyAny=@ 1 ),C[ 1 ](=AnyAny [IS=0,S,N] Any=: Any@ 1 )])(C[NOT=C=0 
,C=S ) C:Any,Any[IS=C ) N]C(=N)N,C[l](Any[IS=0,S,N]Any=:AnyAny=@l 
),C[l](=AnyAny[IS=0,S,N]Any=:Any@l)])(C[NOT=C=0,C=S,C:Any,Any 
[IS=C 5 N]C(=N)N ) C[l](Any[IS=0 ) S,N]Any=:AnyAny-@l),C[l](=AnyAny[ 
IS=0,S,N1 Any=: Any® 1)1) : 


Amines attached to heteroaromatics 


fp2a 


Primary 


C[l](Any[IS=0,S,N]Any=:AnyAny=@l)NH2 
Cril(=AnyAnyriS=0,S,NlAny=:Any@l)NH2 


fp2b 


Secondary 


C[l](Any[IS=0 I S ) N]Any=:AnyAny=@l)NHAny[NOT=H*] 
Cril(=AnyAnyriS=0,S,NlAny=:Any®l)NHAnyfNOT=H*l 


fp2c 


Tertiary 


C[l](Any[IS=0,S,N]Any=:AnyAny=@l)N(Any[NOT=H*])Any[NOT=H*] 
Cril(=AnyAnyriS=0,S,NlAny=:Any@l)N(AnyrNOT=H*])AnyrNOT=H*] 


Anilines 


fp3a 


Primary 


NH2C(:Any)(:Any[NOT=H*l) 


fp3b 


Secondary 


NH(C(:Any)(:AnyrNOT=H*l))Any[NOT=H*l 


fp3c 


Tertiary 


N(C(:Any)(:Any[NOT=H*]))(Any[NOT=H*])Any[NOT=H*] 
N(C(:Any)(:Any[NOT=H*l))=C 
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N in aromatics 




fp4a 


6-membered ring 


Any|"is=N,Cl:N:Anyris=N,C] 




fp4b 


-NH- in heteroaromatics 


N[ 1 ]HAny [IS=C,N1=: Any riS=C,Nl Any[IS=C,N]=: Any[IS=C,N]-@ 1 




fp4c 


-N- in heteroaromatics 


NrilAnyriS=C,Nl=:AnyriS=C,N]AnyriS=C } Nl=:AnyriS=C,Nl-@l 




fp4d 


-N= in heteroaromatics 


N[ 1 ](Any[IS=0,S,N] Any=: AnyAny=@l ) 
N[l](=AnyAny[IS=0,S,N]Any=:Any@l) 




Imines/ Nitrile/ Guanidine/ Amidine 




fp5a 


Imines 


Any [IS=C,H,S]N[NOT=N[ 1 ] (Any [IS=0,S,N] Any=: AnyAny=@ 1 ),N[ 1 ] (=A 
nyAny[IS=0,S,N]Any=:Any@l)]=CfNOT=AnyriS=C,NlC(=N)Nl 




fp5b 


Nitrile 


C#N 




fp5c 


Guanidine 


N[NOT=C[l](Any[IS=0,S,N]Any=:AnyAny=@l)N,C[l](=AnyAny[IS=0,S 
,N]Any=:Any@l)N]C(=N)N[NOT=C[l](Any[IS=0,S,N]Any=:AnyAny=@ 
1 )N,Cr 1 1 (=AnyAny[IS=0,S,N] Any=: Any® 1 )N1 




fp5d 


Amidine (not hetero- 
aromatics) 


Any[NOT=N]C(=N[NOT=N[ 1 ](Any [IS=0,S,N] Any=: AnyAny=@ 1 ),N[ 1 ](= 
Any Any [IS=0,S,N] Any=: Any@ 1 )])N 




N~0/NUro/N=N/N-N 




fp6a 


Hydroxyamine, Oxime, 
Hydroxamic acid....) 


N[!r]-0[!r] 




fp6b 


Nitro, Nitroso 


N(=0) 


m 


fp6c 


N=N Azo (not in a ring) 


N=N[NOT=N[ 1 ](Any pS=0,S,N] Any=: AnyAny=@ 1 ),N[ 1 ](=AnyAny [IS=0 
J S,NlAny=:Any@l)] 


rn 


fp6d 


N-N Hydrazine 


N-N[NOT=N[ 1 ] (Any [IS=0,S,N] Any=: AnyAny=@ 1 ),N[ 1 ] (=AnyAny 
[IS=0,S,N]Any=:Any@l)l 




Amide/ Thioamide/ Sulfonamide 


= 


fp7a 


Amide 1 (NH 2 -CO) 


NH2C=0 




fp7b 


Amide2 (R,-NH-CO) 


Any[NOT=H*]NHC=0 


y * 


fp7c 


Amide3 (R^N-CO) 


AnyrNOT=H*lN(C=0) AnyfNOT=H*] 




fp7d 


Thioamide 1 (NH 2 -CS) 


NH2C-S 


" ! H" 


fp7e 


Thioamide2 (R,-NH-CS) 


AnyfNOT=H*lNHC=S 




fp7f 


Thioamide3 (R,R 2 N-CS) 


Any[NOT=H*]N(C=S)Any[NOT=H*] 


H" ■ 


fp7g 


Sulf.amidel (NH 2 S0 2 ) 


NH2S(=0)(=0) 




fp7h 


Sulf.amide2 (R r NHS0 2 ) 


NH(S(=0)=0)AnyfNOT=H*l 




fo7i 


Sulf.amide3 (R,R r NS0 2 ) 


N(S(=0)=0)(AnyrNOT= H*l)AnyfNOT=H*l 




Alcohol/ Ether/ Aldehyde/ Ketone/ Ester/ Carboxylic acid/ Carbothioic acid/ Sulfinic 




acid/ Sulfonic acid 




fp8a 


Alcohol 


CrNOT=C=0 3 C=S](OH) 




fp8b 


Ether 


Any[NOT=C=0, H*]-0-Any[NOT=C=0,H*l 




fp8c 


Aldehyde 


CCH(=0) 




fp8d 


Ketone 


CC(=0)C 




fp8e 


Ester 


C(=0)OC 




fp8f 


Carboxylic acid 


C(=0)(OH) 




fp8g 


Carbothioic O acid 


C(=S)(OH) 




fp8h 


Carbothioic S acid 


C(=0)(SH) 




fp8i 


sulfinic acid 


Anyris=H,C]S[NOT=S(=0)(=0)] (=0)(OH) 




fp8j 


sulfonic acid 


Any ris=H,Cl S(=0)(=0)(OH) 
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Halogen 


fp9a 


Fluoro 


F 


fp9b 


Chloro 


CI 


fp9c 


Bromo 


Br 


fp9d 


Iodo 


I 


Total C/H/N/ O/ S 


fplOa 


total C 


C 


fplOb 


total H 


H 


fplOc 


total N 


N 


fplOd 


total O 


O 


fplOe 


totals 


S 



